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The seesaw enhancement of the bi-large mixings are discussed for the two- zero textures of 



the neutrino mass matrix. There are no large mixings in both Dirac neutrino mass matrix 
and right-handed Majorana neutrino mass matrix Mr, however, the bi-large mixing is realized 
via the seesaw mechanism. We present twelve sets of and Mr for the seesaw enhancement 
and discuss the related phenomena, the fi — > e + 7 process and the leptogenesis. The decay 
rate of fi — > e + 7 is enough suppressed due to zeros in the Dirac neutrino mass matrix. Six mo 
lead to the lepton asymmetry, which can explain the baryon number in the universe. Other 
six tud are the real matrices, which give no CP asymmetry. Modified Dirac neutrino mass 
matrices are also discussed. 
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The texture with zeros of the neutrino mass matrix have been discussed [1, 2, 3, 4] to 
explain neutrino masses and mixings [5], which have been presented by the recent neutrino 
experiments [6, 7, 8, 9]. It was found that the two- zero textures are consistent with the 
experimental data in the basis of the diagonal charged lepton mass matrix [10]. Consequently, 
the neutrino mass matrix does not display the hierarchical structure as seen in the quark mass 
matrix [11, 12, 13, 14, 15, 16, 17]. 

Since the two-zero textures of ref . [10] are given for the light effective neutrino mass matrix 
M„, one needs to find the seesaw realization [18] of these textures from the standpoint of the 
model building. We have examined the seesaw realization of the neutrino mass matrix with two 
zeros [19]. Without fine tunings between parameters of the Dirac neutrino mass matrix mp 
and the right-handed Majorana neutrino one M R , we obtained several textures of m D for the 
fixed M R [19]. Among them, there are textures of m D and M R which have hierarchical masses 
without large mixings. These present the seesaw enhancement of mixings, because there is no 
large mixings in mo and Mr, but the bi-large mixing is realized via the seesaw mechanism. 
The seesaw enhancement are important in the standpoint of the quark-lepton unification, in 
which quark masses are hierarchical and quark mixings are very small 1 . 

The general discussions of the seesaw enhancement were given in the case of two flavors 
[20, 21]. Specific cases were discussed in the case of three flavors [22, 23] because it is very 
difficult to get general conditions for the seesaw enhancement of the bi-large mixing. 

However, the two-zero texture of the neutrino mass matrix M v are helpful to study the 
seesaw enhancement of the bi-large mixing. In this paper, we present sets of mc and Mr to 
give the seesaw enhancement in the two-zero textures of M v and discuss the related phenomena, 
the (/, — > e + 7 process and the leptogenesis [24] . 

There are fifteen two-zero textures for the neutrino mass matrix M v , which have five inde- 
pendent parameters. Among these textures, seven acceptable textures with two independent 
zeros were found for the neutrino mass matrix [10], and they have been studied in detail 

[12, 13, 16, 17]. Especially, the textures A 1 and A 2 of ref. [10], which correspond to the hierar- 

1 Although phcnomcnological analyses of the two-zero textures were given in the diagonal basis of the charged 
lepton, some authors [15, 16, 17] have also studied the two-zero textures of neutrinos in the basis of charged 
lepton mass matrix with small off-diagonal components. 
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chical neutrino mass spectrum, are strongly favored by the recent phenomenological analyses 
[11, 12, 17]. Therefore, the two textures are taken in order to discuss the seesaw enhancement. 
Putting data of neutrino masses and mixings [25] , 

0.35 < tan 2 9 sun < 0.54 , 6.1 x 10~ 5 < Am s 2 un < 8.3 x 10^ 5 eV 2 , 90%C.L. , 

0.90 < sin 2 26 atm , 1.3 x 10~ 3 < Am 2 tm < 3.0 x 10~ 3 eV 2 , 90%C.L. , (1) 

the relative magnitude of each entry of the neutrino mass matrix is roughly given for the 

textures A 1 and A 2 as follows: 

/0 A\ /OA 0\ 

M v ~ m 1 1 for A x , m A 1 1 for A 2 , (2) 

\\ 1 l) Vo 1 \) 

where m denotes a constant mass and A ~ 0.2. These matrix is given in terms of the Dirac 
neutrino mass matrix m D and the right-handed Majorana neutrino mass matrix M R by the 
the seesaw mechanism as 

M v = m D M R X m T D . (3) 

Zeros in mo and M R provide zeros in the neutrino mass matrix M v of eq.(2) as far as we 
exclude the possibility that zeros are originated from accidental cancellations among matrix 
elements. In other words, we take a standpoint that the two-zero texture should come from 
zeros of the Dirac neutrino mass matrix and the right-handed Majorana mass matrix. Possible 
textures of m D and M R were given in ref. [19]. Among them, we select the set of and M R , 
which reproduce the seesaw enhancement of the bi-large mixing. 

Let us fix the right-handed Majorana neutrino mass matrix without large mixings. We take 
simple right-handed Majorana neutrino mass matrix with only three independent parameters. 
Then, there are = 20 textures. Among them, six textures are excluded because they have 
a zero eigenvalue, which corresponds to a massless right-handed Majorana neutrino. Other 
two textures are also excluded because the two-zero textures A\ and A 2 cannot be reproduced 
without accidental cancellations. One of the two textures is the diagonal matrix, and another 
one is the matrix with three zeros in the diagonal elements. 

We show twelve real mass matrices with three independent parameters 2 with mass eigen- 
values I Mi I = A m M 3 and \M 2 \ = A n M 3 , where M 3 is the mass of the third generation, and m 
2 The classification of Mr, types a^, bi, a (i — 1,2,3,4), follows from ref.[19]. 



and n are integers with m > n > 1: 
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where there are no large mixings in twelve matrices since the mass eigenvalues are supposed 
to be hierarchical. The minus signs in the matrix elements are taken to reproduce signs in the 
texture A 1 and A 2 of eq.(2). 

There are several Dirac neutrino mass matrices to give the textures A\ and A 2 in eq.(2) 
[19]. We show Dirac neutrino mass matrices (mo) a*, ( m D)bi, (jno)^ (2 = 1 ~ 4) with maximal 
number of zeros, which have no large mixings, to give the texture A 2 3 . For each matrix of 
(M R ) a ., (M R ) bi , (M R ) C . those are given as follows: 

/A 

ai type m D ~ m D0 





m D0 
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3 For the texture Ai, we easily obtain the Dirac neutrino mass matrices by exchanging the second and third 
rows. 
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bi type m D ~ 
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where m^o denotes the magnitude of the Dirac neutrino mass and complex coefficients of 
order one are omitted. Although these matrices have no large mixing among three families, 
the neutrino mass matrix M v has the bi-large mixing through the seesaw mechanism. These 
are so called seesaw enhancement of the bi-large mixing. 

These Dirac matrices are ones with maximal number of zeros. Without changing the 
mixings and the mass eigenvalues in the leading order, some zeros can be replaced with small 
non-zero entries as follows: 
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Ci type m D ~ m D0 \ \ x \ y A 2 ] , > -j 
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where x, y and z are positive integers. These Dirac neutrino mass matrices are asymmetric 
ones. However, only the 63 and 64 textures in eq.(ll) are adapted to the symmetric texture in 
the S'O(10)-like GUT if y, m and n are relevantly chosen [15, 14]. For example, in the 63 case, 
taking y = n/2 and m = n + 2, the symmetric mass matrix is given, especially, putting n — 8, 
we have the hierarchical mass matrix such like the up-quark mass matrix. 

Let us discuss these obtained textures of mp in the \i — > e + 7 decay and the leptogenesis. It 
is well known that the Yukawa coupling of the neutrino contributes to the lepton flavor violation 
(LFV). Many authors have studied the LFV in the minimal supersymmetric standard model 
(MSSM) with right-handed neutrinos assuming the relevant neutrino mass matrix [26, 27, 28, 
29, 30, 31]. In the MSSM with soft breaking terms, there exist lepton flavor violating terms 
such as off-diagonal elements of slepton mass matrices and trilinear couplings (A-term). It 
is noticed that large neutrino Yukawa couplings and large lepton mixings generate the large 
LFV in the left-handed slepton masses. For example, the decay rate of /i — > e + 7 can be 
approximated as follows: 



e 2 

e 5 



r(/i^e + 7 ) ~ j^mlF 



(6 + 2a^| (y ^ f)2iln 



M 



x 



IGtt 2 11 ^ M r ' (13) 

where the neutrino Yukawa coupling matrix Y v is given as Y v = m D /v 2 (v 2 is a VEV of Higgs) at 
the right-handed mass scale Mr, and F is a function of masses and mixings for SUSY particles. 
In eq.(13), we assume the universal scalar mass (m 50 ) for all scalars and the universal A-term 
(Af — a m S oYf) at the GUT scale M X - Therefore the branching ratio /x — > e + 7 depends 
considerably on the texture m D [29, 30, 31]. 

The magnitude of (ra 0^0)21 is a key ingredient to predict the branching ratio of the 
/i — > e + 7 process 4 . Many works have shown that this branching ratio is too large [29, 30]. 
4 mDm^ does not depend on the basis of the right-handed sector. 
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The conditions for (rriDm D ) 21 were given in ref. [32] as follow: 
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where 



#ij = Y,( m D)ik{m ] D ) kj \n 
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(15) 
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These conditions give constraints for the magnitude of [m D rrv D )ij and M 3 . Zeros in Dirac 
mass matrices m D may lead to (m D m} D )ij = and then it suppress the \i — * e + 7 decay. 
Actually, all mo in eqs.(7), (8) and (9) give (morn 0)21 = and (m^m^, ) 31 (m^m^,) 23 = 5 . 
Even if non-zero terms \ x , \ y , \ z are taken as seen in eqs.(10), (11) and (12), (m^m^^i and 
{m D m} D ) 31 (m D m^ D ) 23 are suppressed as far as x,y,z ^> 1. Then, the branching ratio is safely 
predicted to be below the present experimental upper bound 1.2 x 1CT 11 [33] due to zeros. 

Let us examine our textures in the leptogenesis [34, 35, 36], which is based on the Fukugita- 
Yanagida mechanism [24]. The CP violating phases in the Dirac neutrino mass matrix are key 
ingredients for the leptogenesis while the right-handed Majorana neutrino mass matrix are 
taken to be real in eqs.(4), (5) and (6). Although the non-zero entries in the Dirac neutrino 
mass matrix are complex, three phases are removed by the re-definition of the left-handed 
neutrino fields. There is no freedom of re-definition for the right-handed ones in the basis with 
the real Mr. We should move to the diagonal basis of the right-handed Majorana neutrino 
mass matrix in order to calculate the magnitude of the leptogenesis. Then, the Dirac neutrino 
mass matrices mo in the new basis is given as follows: 



where Pl is a diagonal phase matrix and Or is the orthogonal matrix which diagonalizes Mr 
as OrMrOr in eqs.(4), (5) and (6). Since the phase matrix Pl can remove one phase in each 
row of rriD, three phases disappear in m D . 



m D = P L m D Or , 



(16) 



5 For the texture A\ case, the some Dirac mass matrices give non-zero (m£>ra£>)2i, which leads to the 
constraint for M3 
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As a typical example, we show the case of the 63 texture in eq.(5). By taking three eigen- 
values of M R as follows 6 : 

M 1 = X m M 3 , M 2 = -X n M 3 . (17) 

We obtain the orthogonal matrix Or as 

(cos 9 sin # \ 
-sin# cos 9 , tan 2 9 = \ m ~ n . (18) 

l) 

Then the Dirac mass matrices of 63 in eq.(8) can be parameterized in the new basis as 
follows: 



/ Af +1 0\ 



m D = m D0 



Or , (19) 



At 
\ \*e ip I J 

where only one phase p remains. The magnitude of moo is determined by the relation m 2 m 



m M 3 , where m ~ ^Aml tm /2. 

We examine the lepton number asymmetry in the minimal SUSY model with the right- 
handed neutrinos. In the limit Mi <C M 2 , M 3 , the lepton number asymmetry ei (CP asym- 
metry) for the lightest heavy Majorana neutrino (Ni) decays into / T ± [37] is given by 

where v 2 = v sin (3 with v = 174GeV. The lepton asymmetry Y L is related to the CP asymmetry 
through the relation 

n = 2£Z»E = , (21) 
s g* 

where s denotes the entropy density g* is the effective number of relativistic degrees of freedom 
contributing to the entropy and k is the so-called dilution factor which accounts for the washout 
processes (inverse decay and lepton number violating scattering). In the MSSM with right- 
handed neutrinos, one gets g* = 232.5. 

The produced lepton asymmetry Yl is converted into a net baryon asymmetry Yb through 
the (B + L)-violating sphaleron processes. One finds the relation [38] 

£ 8 N f + 4 N H 



6 The minus sign of Mi is necessary to reproduce Mr in eq.(5). This minus sign is transferee! to mo by the 
right-handed diagonal phase matrix diag(l,i, 1). 



where Nf and Nh are the number of fermion families and Higgs doublets, respectively. Taking 
into account Nf = 3 and N H = 2 in the MSSM, we get 

Y B = ~Y L . (23) 

On the other hand, the low energy CP violation, which is a measurable quantity in the 
long baseline neutrino oscillations [39], is given by the Jarlskog determinant J C p [40], which is 
calculated by 

det[M e Mj,M u Ml] = -2iJ C p(ml-ml)(ml-ml)(ml-ml)(ml-ml)(ml-mj)(ml-ml), (24) 

where Mn is the diagonal charged lepton mass matrix, and mi, m 2 , m 3 are neutrino masses. 

Since the CP violating phase is only p, we can find a link between the leptogenesis (e x ) and 
the low energy CP violation (Jcp) in our textures of the Dirac neutrinos. By using the Dirac 
neutrino mass matrix in eq.(19), we get 

3777 2 

ei ~ ^ ^ m sin2p ~ -8.8 x 10~ 17 Mi sin2p , 

8ttt;2 

where Mi is given in the GeV unit and tan/3 > 10 is taken. It is remarked that e x only 
depends on Mi and the phase p, and the relative sign of e x and Jcp is opposite. Taking the 
experimental data Am s 2 ol / Am 2 tm ~ A 2 and sin2p ~ 1, we predict Jcp — 0.01, which is rather 
large and then is favored for the future experimental measurement. 

The five cases of the Dirac neutrino mass matrix (ai, a 3 , bi, c\, c 3 ) in eqs.(7), (8) and 
(9) lead to same results in eq.(25). In other six cases of the Dirac neutrino mass matrix 
(02,04,62,64,02,04), the CP violating phases are removed because of only three non-zero en- 
tries. Then, we get ei = 0, but the same result in eq.(25) for Jcp- 

If we use the modified Dirac neutrino mass matrices in eqs.(10), (11) and (12), new CP 
violating phases appear. However, the contribution to e x is a next-leading one as far as 1 > 
1, y > l,z > 1. 

In order to calculate the baryon asymmetry, we need the dilution factor involves the inte- 
gration of the full set of Boltzmann equations [41]. A simple approximated solution which has 
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k — 0.3 ( ( in ) . (iO" a eV < mi < 10 3 eV) , (26) 



been frequently used is given by [42] 

- 3 eV \ / mi ^ - 6 
mi J V 11 10" 3 eV 

where 

mi = Mi • (27) 

By using this approximate dilution factor and eqs.(21) and (22), we can estimate Yb in our 
textures as follows: 

Y B ~ -2.3 x ei re . (28) 

It is noticed that Yb and Jcp are same sign since e\ has minus sign. 
The WM AP has given the new result [43] 

VB = 6.5+^ x 10- 10 (1 a) , (29) 

which leads to 

Y B ^ \vb ■ (30) 

In our textures, we have (m.D*mD)ii = m 2 D0 X m , which gives mi = \ \J Aml tm ~ 0.022. Then 
we get the dilution factor k~7x 10~ 3 . Putting the observed value into eq.(28), we get 

M 1 sin 2p ~ 6 x 10 10 GeV . (31) 

This result means that Mi is should be larger than 6 x 10 10 GeV in order to explain the baryon 
number in the universe. This value is consistent with previous works [34, 35, 36]. 

It is important to present the discussion from the standpoint of the GUT, which is given 
after eq.(ll). Taking n = 8 and m = 6 in the 63 case of eq.(ll) as in the previous discussion, 
one obtains M 3 ~ 10 15 GeV and M x ~ 10 8 GeV taking account of Am^ tm ~ 2 x 10~ 3 eV 2 . This 
result does not satisfy the condition of eq.(31). However, the simple 50(10) fermion mass 
relation may be consistent with the leptogenesis in the case of the more complicated texture of 
Mr, which leads to the two- zero texture A 2) as seen in the work of [44]. Details are presented 
in the preparing paper including the degenerate case of Mr in the simple 50(10) approach 
[45]. 
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We add the discussion of another important problem. In the framework of supersymmetric 
thermal leptogenesis, there is cosmological gravitino problems. The gravitino with a few TeV 
mass does not favor Mi > 10 10 GeV [46], because Mi should be lower than the maximum 
reheating temperature of the universe after inflation. In order to keep the thermal leptogenesis 
in the SUSY model, we may consider the gravitino with O(100)TeV mass, which is derived 
from the anomaly mediated SUSY breaking mechanism [47]. 

Summary is given as follows. We have discussed the textures with the seesaw enhancement. 
These textures are important in the standpoint of the quark-lepton unification, in which quark 
masses are hierarchical and quark mixings are very small. It is very difficult to get general 
conditions for the seesaw enhancement of the bi-large mixing, however, the two- zero texture of 
the left-handed neutrino mass matrix M„ are helpful to study the seesaw enhancement of the 
bi-large mixing. Once the basis of the right-handed Majorana neutrino mass matrix is fixed, 
one can find some sets oirriD and Mr, which have hierarchical masses without large mixings, to 
give the two- zero textures A 1 and A 2 without fine tuning among parameters of these matrices. 
These sets present the seesaw enhancement of the bi-large mixing, because there is no large 
mixings in and Mr, but bi-large mixing is realized via the seesaw mechanism. We present 
twelve sets of m D and Mr for the seesaw enhancement. Then, the decay rate of ji — ■> e + 7 
is enough suppressed due to zeros in the Dirac neutrino mass matrix. Six sets lead to the 
lepton asymmetry, which depends on only Mi and the phase p. Putting the observed value 
of baryon number in the universe, M x ~ 6 x 10 10 GeV is obtained. It is remarked that J C p is 
the same sign as the Y B , and its magnitude is predicted to be ~ 0.01. Other six ones provide 
the real Dirac neutrino mass matrices, which give no CP asymmetry. Study of modified right- 
handed Majorana neutrino mass matrices is important for realistic model buildings based on 
the quark-lepton unification. 

This research was supported in part by the National Science Foundation in the Grant-in-Aid 
for Science Research, Ministry of Education, Science and Culture, Japan(No. 12047220). 



11 



References 

[I] H. Nishiura, K. Matsuda and T. Fukuyama, Phys. Rev. D60 (1999) 013006. 

[2] E. K. Akhmedov, G. C. Branco, M. N. Rebelo, Phys. Rev. Lett. 84 (2000) 3535. 

[3] S.K. Kang and C.S. Kim, Phys. Rev. D63 (2001) 113010. 

[4] K.T. Mahanthappa and Mu-C. Chen, Phys. Rev. D62 (2000) 113007. 

[5] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28 (1962) 870. 

[6] Super-Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 81 (1998) 1562; ibid. 
82 (1999) 2644; ibid. 82 (1999) 5194. 

[7] Super-Kamiokande Collaboration, S. Fukuda et al. Phys. Rev. Lett. 86 (2001) 5651; 5656. 

[8] SNO Collaboration: Q. R. Ahmad et al., Phys. Rev. Lett. 87 (2001) 071301; 
Phys. Rev. Lett. 89 (2002) 011301; ibid. 89 (2002) 011302; nucl-ex/0309004. 

[9] KamLAND Collaboration, K. Eguchi et al., Phys. Rev. Lett. 90 (2003) 0212021. 

[10] PH. Frampton, S.L. Glashow and D. Marfatia, Phys. Lett. B536 (2002) 79. 

[II] Z. Xing, Phys. Lett. B530 (2002) 159. 

[12] W. Guo and Z. Xing, Phys. Rev. D67 (2003) 053002. 

[13] R. Barbieri, T. Hambye and A. Romanino, JHEP 0303 (2003) 017. 

[14] K.T. Mahanthappa and Mu-C. Chen, Phys. Rev. D65 (2002) 053010; D68 (2003) 017301. 

[15] M. Bando and M. Obara, Prog. Theor. Phys. 109 (2003) 995; hep-ph/0212242; 
M. Bando, S. Kaneko, M. Obara and M. Tanimoto, hep-ph/0309310. 

[16] A. Ibara and G. G. Ross, Phys. Lett. B575 (2003) 279. 

[17] M. Honda, S. Kaneko and M. Tanimoto, JHEP 0309 (2003) 028. 



12 



[18] Gell-Mann, P. Ramond and R. Slansky, in Supergravity, Proceedings of the Workshop, 
Stony Brook, New York, 1979, edited by P. van Nieuwenhuizen and D. Freedmann, North- 
Holland, Amsterdam, 1979, p. 315; 

T. Yanagida, in Proceedings of the Workshop on the Unified Theories and Baryon Number 
in the Universe, Tsukuba, Japan, 1979, edited by O. Sawada and A. Sugamoto, KEK 
Report No. 79-18, Tsukuba, 1979, p.95. 

[19] A. Kageyama, S. Kaneko, N. Shimoyama and M. Tanimoto, Phys. Lett. B538 (2002) 96. 

[20] A. Yu. Smirnov, Phys. Rev. D48 (1993) 3264. 

[21] M. Tanimoto, Phys. Lett. 345B (1995) 477. 

[22] M. Bando, T. Kugo and K. Yoshioka, Phys. Rev. Lett. 80 (1998) 3004; 

M. Jezabek and Y. Sumino, Phys. Lett. B440 (1998) 327; ibid. B457 (1999) 139. 

[23] G. Altarelli, F. Ferruglio and I. Masina, Phys. Lett. B472 (2000) 382; 

E.K. Akhmedov, G.C. Branco and M.N. Rebelo, Phys. Lett. B478 (2000)215; 
A. Datta, F.-S. Ling, and P. Ramond, Nucl. Phys.B671 (2003) 383; 
S. Lavignac, I. Masina and C. A. Savoy, Nucl. Phys.B633 (2002) 139; 
W. Rodejohann, hep-ph/03 11142. 

[24] M. Fukugita and T. Yanagida. Phys. Lett. B175 (1986) 45. 

[25] G. L. Fogli, E. Lisi, M. Marrone, D. Montanino, A. Palazzo and A.M. Rotunno, Phys. 
Rev. D67 (2003) 073002; 

J. N. Bahcall, M. C. Gonzalez-Garcia and C. Peha-Garay, JHEP 0302 (2003) 009; 

M. Maltoni, T. Schwetz and J.W.F. Valle, Phys. Rev. D67 (2003) 093003; 

PC. Holanda and A. Yu. Smirnov, JCAP 0302 (2003) 001; 

V. Barger and D. Marfatia, Phys. Lett. 555B (2003) 144; 

M. Maltoni, T. Schwetz, M. Tortola and J.W.F. Valle, hep-ph/0309130. 

[26] F. Borzumati and A. Masiero, Phys. Rev. Lett. 57 (1986) 961. 



13 



[27] J. Hisano, T. Moroi, K. Tobe, M. Yamaguchi and T. Yanagida, Phys. Lett. B357 (1995) 
579; 

J. Hisano, T. Moroi, K. Tobe and M. Yamaguchi, Phys. Rev. D53 (1996) 2442. 

[28] J. Hisano, D. Nomura and T. Yanagida, Phys. Lett. B437 (1998) 351; 
J. Hisano and D. Nomura, Phys. Rev. D59 (1999) 116005; 

M.E. Gomez, G.K. Leontaris, S. Lola and J. D. Vergados, Phys. Rev. D59 (1999) 116009; 

W. Buchmiiller, D. Delepine and F. Vissani, Phys. Lett. B459 (1999) 171; 

W. Buchmiiller, D. Delepine and L.T. Handoko, Nucl. Phys. B576 (2000) 445; 

J. Ellis, M.E. Gomez, G.K. Leontaris, S. Lola and D.V. Nanopoulos, Eur. Phys. J. C14 

(2000) 319; 

J. L. Feng, Y. Nir and Y. Shadmi, Phys. Rev. D61 (2000) 113005; 

S. Baek, T. Goto, Y. Okada and K. Okumura, Phys. Rev. D63 (2001) 051701. 

[29] J. Sato, K. Tobe, and T. Yanagida, Phys. Lett. B498 (2001) 189; 
J. Sato and K. Tobe, Phys. Rev. D63 (2001) 116010; 
S. Lavignac, I. Masina and C. A. Savoy, Phys. Lett. B520 (2001) 269. 

[30] J. A. Casas and A. Ibarra, Nucl. Phys. B618 (2001) 171. 

[31] A. Kageyama, S. Kaneko, N. Simoyama and M. Tanimoto, Phys. Lett. B527 (2002) 206; 
Phys. Rev. D65 (2002) 096010. 

[32] J. Ellis, J. Hisano, M. Raidal and Y. Shimizu, Phys. Rev. D66 (2002) 115013. 

[33] MEGA Collaboration, M. L. Brooks et al., Phys. Rev. Lett. 83 (1999) 1521. 

[34] M. Flanz, E. A. Paschos and U. Sarkar, Phys. Lett. B345 (1995) 248; 

W. Buchmiiller and M. Pliimacher, Phys. Lett. B389 (1996) 73; ibid. B431 (1998) 354; 

Int. J. Mod. Phys. A15 (2000) 5047; 

A. Pilaftsis, Phys. Rev. D56 (1997) 5431. 

[35] A. S. Joshipura, E. A. Paschos and W. Rodejohann, JHEP 0108 (2001) 029; 
W. Buchmiiller and D. Wyler, Phys. Lett. B 521 (2001) 291; 

14 



J. Ellis and M. Raidal, Nucl. Phys. B643 (2002) 229; 
S. Davidson and A. Ibarra, Nucl. Phys. B648 (2003) 345; 

G. Branco, T. Morozumi, B. Nobre and M.N. Rebelo, Nucl. Phys. B 617 (2001) 475; 

G. Branco, R. Gonzalez Felipe, F. R. Joaquim and M. N. Rebelo, Nucl. Phys. B 640 (2002) 
202; 

P. Frampton, S. Glashow and T. Yanagida, Phys. Lett. B548 (2002) 119; 

T. Endoh, S. Kaneko, S. K. Kang, T. Morozumi and M. Tanimoto, Phys. Rev. Lett. 89 

(2002) 231601; 

S. Kaneko, M. Katsumata and M. Tanimoto, JHEP 0307 (2003) 025. 

[36] D. Falcone, Phys. Rev. D66 (2002) 053001; 
Z. Xing, Phys. Lett. B545 (2002) 352; 
W. Rodejohann, Phys. Lett. B542 (2002) 100; 

M. Fujii, K. Hamaguti and T. Yanagida, Phys. Rev. D65 (2002) 115012. 

[37] L. Covi, E. Roulet and F. Vissani, Phys. Lett. B 384 (1996) 169. 

[38] J. A. Harvey and M. S. Turner, Phys. Rev. D42 (1990) 3344; 

S. Y. Khlebnikov and S. E. Shaposhnikov, Nucl. Phys. B308 (1998) 169. 

[39] M. Tanimoto, Phys. Rev. D55 (1997) 322; Prog. Theor. Phys. 97 (1997) 901; Phys. Lett. 
B435 (1998)373; Phys. Lett. B462 (1999) 115; 
J. Arafune and J. Sato, Phys. Rev. D55 (1997)1653; 
J. Arafune, M. Koike and J. Sato, Phys. Rev. D56 (1997) 3093; 

H. Minakata and H. Nunokawa, Phys. Lett. B413 (1997) 369; Phys. Rev. D57 (1998) 4403. 
For recent proposal of the search for CP violation in neutrino oscillation, see Y. Itow et.al., 
hep-ex/0106019. 

[40] C. Jarlskog, Phys. Rev. Lett. 55 (1985) 1039. 

[41] M. Luty, Phys. Rev. D45 (1992) 455; 
M. Plumacher, Z. Phys. C74 (1997) 549. 



15 



[42] E. W. Kolb, M. S. Turner, The early universe, Redwood City, USA: Addison- Wesley 
(1990), (Frontiers in physics, 69); 

M. Flanz and E. A. Paschos, Phys. Rev. D58 (1998) 113009; 

A. Pilaftsis, Int. J. Mod. Phys. A14 (1999) 1811; 

H. B. Nielsen and Y. Takanishi, Phys. Lett. B 507 (2001) 241. 

[43] C. L. Bennett et al, Astrophys.J. Suppl. 148 (2003) 1; 
D. N. Spergel et al., Astrophys.J. Suppl. 148 (2003) 175; 
H. V. Peiris et al, Astrophys.J. Suppl. 148 (2003) 213. 

[44] W. Buchmiiller and M. Pliimacher, Phys. Lett. B521 (2001) 291. 

[45] M. Bando, S. Kaneko, M. Obara and M. Tanimoto, to be appeared in hep-ph. 

[46] M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D63 (2001) 103502; astro-ph/0402490. 

[47] L. Randall and R. Sundrum, Nucl. Phys. B557 (1999) 79; 

G.F. Giudice, M. A. Luty, H. Murayama and R. Rattazzi, JHEP 9812 (1998) 027. 



16 



